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ON THE DYNAMICS OF SHORT PRESSURE PROBES: SOME DESIGN 


FACTORS AFFECTING FREQUENCY RESPONSE 
by Ted W. Nyland, David R. Englund, and Robert C. Anderson 
Lewis Research Center 

SUMMARY 

Frequency response measurements of 2. 54 -centimeter -long probes with 0. 325- and 
0. 160 -centimeter diameters are described. Comparison of the measured response and 
that predicted from existing analyses showed good correlation. The results indicate that 
the Bergh and Tijdeman equations accurately predict the response of probes which in- 
clude extension tubes for measuring time -average pressure. It is shown that the use of 
some type of baffle for particle impact protection which involves an enlargement of the 
entrance flow area does not reduce the resonant frequency of the probe. A computer pro- 
gram for calculation of probe response based on the Bergh and Tijdeman equation is 
given in the appendix. 


INTRODUCTION 

Investigations of the dynamic stall characteristic of turbojet engines have placed 
great emphasis on the measurement of transient total pressure. Experimenters wish to 
determine how and at what rate stall propagates through an engine. To obtain this in- 
formation, many transient total pressure probes are located in the inlet and in between 
stages of the compressor on the engine under test. In general, the requirements placed 
on these transient measuring systems include a frequency -response flat (amplitude ratio 
within ±5 percent) to at least 500 -hertz, operation in airstream temperatures ranging 
from room temperature to approximately 800 K and reliable operation where particulate 
matter may be entrained in the airstream. In addition, it is of interest to measure the 
average total pressure at each probe location. The average pressure may range be- 
tween a few to 200 newtons per square centimeter, and the transient pressure amplitudes 
may reach 30 percent of these values. 

Measuring systems used in these investigations have consisted of probes with inter- 



nally mounted miniature pressure transducers connected to short total head tubes (ref. 1). 
Such probes are generally water cooled so that the transducer is maintained at a mod- 
erate and relatively constant temperature. To protect the transducer from damage by 
the entrained particles in the airstream, the transducers have been mounted at right 
angles to the axis of the probe, or offsets or baffles have been included in the internal 
geometry of the probe. The measurement of average total pressure has been accom- 
plished by connecting each probe through small -diameter tubing to high -accuracy low- 
frequency -response pressure instrumentation located external to the engine. This 
approach has proved to be preferable to relying on the average value of the output of the 
miniature (transient) pressure transducers mainly because the zero stability of these 
transducers is not compatible with the accuracy required for these data. The problem 
faced by the probe designer is how to incorporate these features of particle impact pro- 
tection and average total pressure measurement into the probe design without degrading 
the transient response of the measuring system. 

The purpose of this report is to present practical design information related to op- 
timizing the transient response of such probes. The information contained herein is the 
result of an experimental program at the Lewis Research Center. The major objectives 
of this program were 

(1) To measure the frequency response of simulated transient pressure probes 

(2) To measure the effect on response of variations in probe geometry aimed at pro- 

viding particle impact protection and average pressure measuring capability 

(3) To correlate measured response with existing theoretical analyses on transient 

response of probes. 

There are a number of theoretical analyses that can be used to predict the transient 
response of simple probe configurations. Estimates of the natural frequency can be ob- 
tained by using the simple quarter wavelength (organ pipe) and other tube -and -volume 
analyses found in many textbooks (e. g. , see refs. 2 and 3). These analyses are based 
on linear second -order models of the system. More rigorous analyses have been done 
by Iberall (ref. 4) and more recently by Bergh and Tijdeman (ref. 5). The Iberall 
analysis is limited to a single tube and volume; Bergh and Tijdeman extended this analy- 
sis to a series connection of a number of tubes and volumes. These later analyses are 
not based on second-order system models and consequently do not yield a simple expres- 
sion for a resonant frequency. However, they do permit the calculation of the amplitude 
ratio and phase shift as a function of excitation frequency from which the frequency for 
maximum amplitude ratio, that is, the resonant frequency, can be obtained. In either 
case the useful frequency range of the probe is that portion of the frequency spectrum 
for which the amplitude ratio is equal to one to within some acceptable tolerance (±5 per- 
cent for our purposes). In cases where the damping is negligible, which is generally 
true for the probes considered in this report, the frequency at the 5-percent tolerance 
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limit is approximately equal to 20 percent of the resonant frequency. 

The analyses of references 4 and 5 have been verified experimentally for frequencies 
up to about 2000 hertz (refs. 5 and 6). Recently, Goldschmied (ref. 7) has further veri- 
fied the analysis of reference 4. The data contained herein show that references 3 and 
4 accurately predict the response of the short tubes tested in this work. 


PROBE CONFIGURATIONS 

The probe designs used in this work were intended to simulate only the internal geo- 
metry of typical total pressure probes. The outside diameter, wall thickness, and flat 
end of the probes were chosen for ease of mounting in the test apparatus and do not rep- 
resent good design practice for total pressure probes. The exterior and inlet geometries 
of probes for actual measurements in a flowing stream should be governed by such con- 
siderations as flow angle variations, pressure gradients, probe blockage, etc. 


Straight probe 

The basic (straight) probe tested in this work consisted of a short straight tube of 
constant inside diameter terminated with a cavity formed by the transducer mounting 
(see fig. 1). As the volume ratio, the ratio of cavity volume to tube volume, is a pa- 
rameter in the analytical equations, provision was made for varying the cavity volume. 

In this report, the straight probes will be designated as either A1 or A2, referring to 
the diameters 0. 325 or 0. 160 centimeter, respectively, followed by a number corres- 
ponding to the volume ratio. A miniature quartz piezoelectric pressure transducer was 
used in this probe. 

In order to provide a means to measure the average total pressure, the probe was 
provided with taps drilled at right angles to the probe axis as shown in figure 1. The tap 
diameter was either 0. 150 or 0. 071 centimeter. These matched the inside diameters of 
tubing used when teeing the effects of average pressure measuring systems on probe 
response. The lengths of the average pressure tubing tested were 15, 61, and 183 centi- 
meters. 

In one variation of the basic probe design, the transducer was mounted at right 
angles to the probe axis as shown in figure 2 . The cavity volume was cylindrically 
shaped with its axis perpendicular to the tube axis. A miniature semiconductor strain 
gage pressure transducer, 0. 317 centimeter in diameter, was used in this probe con- 
figuration. 
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Baffled Probe 


The baffled probe tested in this work was designed to provide particle impact pro- 
tection for the transducer. The baffle was built into the opening of the probe as shown 
in figure 3. The baffle consisted of an annular passageway connected through a set of 
holes to the inside tube. Four similar configurations were tested, each having different 
ratios of annular area to tube area and hole area to tube area as shown in the figure. 

For baffle probes Bl, B2, and B3, a set of 5 holes, 0. 139 centimeter in diameter, 
were drilled between the annulus and probe tube. For probe B4, two sets of holes were 
used. The volume ratio based on inside tube length was 0. 118. The pressure transducer 
was of'the miniature quartz variety. 


TEST APPARATUS AND PROCEDURE 

A sinusoidal pressure generator developed at the Lewis Research Center was used 
to determine the frequency response of the probes. A drawing of the generator is shown 
in figure 4, and a detailed description of the generator is given in reference 8. The gen- 
erator consists of a 2. 54 -centimeter -diameter closed tube driven by an annular shape 
air jet. The frequency of oscillation is varied by moving a tuning piston in the tube, 
effectively changing the tube length. Frequencies between 300 and 5000 hertz can be ob- 
tained with the generator. The probes are flush mounted to the inside surface of the 
resonator tube. A second pressure transducer is flush mounted directly opposite the 
probes and is used as a reference for both the amplitude and phase angle measurements. 
All probes were tested at atmospheric pressure and at room temperature with nominal 
pressure amplitudes of 0. 5 newton per square centimeter (5 percent of atmospheric) 
peak to peak. The estimated error in the amplitude ratios, which are less than 15 deci- 
bels is 2 to 3 percent. As the amplitude ratios increase beyond 15 decibels, the error 
in the ratio increases to 10 to 20 percent at 30 decibels. Estimated error for the fre- 
quency measurement is ±1 hertz and for the phase angle measurement is ±5°. 


TEST RESULTS 

Effects of Volume Ratio on Probe Response 

The measured frequency response of probe A1 -0.027 is shown in figure 5. Also 
shown is the analytical response curve for the probe as calculated from reference 4. 
(The values of fluid parameters used in all calculations are given in table I. ) The ana- 
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lytical curve includes an end correction of 8D/3 jt, as suggested in reference 3, added to 
the length term (D is tube diameter). Without the end correction, the analytical curve 
is the same shape, but it has a resonant frequency approximately 10 percent higher than 
the curve shown. Generally, the measured resonant frequencies of all the A1 and A2 
probes were within 5 percent of the end corrected analytical values. 

The effects of changes in volume ratio on the resonant frequency are shown in fig- 
ure 6. Plotted is the percent reduction in resonant frequency as a function of volume 
ratio as calculated from reference 4 with end correction. The reduction in resonant 
frequency is calculated as the difference between the resonant frequency of a probe with 
(Fy) and without (Fq) a volume at the transducer end. Calculations indicate that the re- 
sonant frequency is slightly dependent on probe diameter. For the two probe diameters 
tested in this program, the resonant frequency for the same length and volume ratio dif- 
fered by approximately 3 percent. For example, Fq for the A1 probe was 3000 hertz 
and for the A2 probe, 3100 hertz. As this is not a significant difference, the analytical 
curve plotted in figure 6 is the average value for the A1 and A2 probes. Also plotted in 
the figure are experimental data from a number of straight probe tests. For volume 
ratios less than 0. 2, the experimental data are within 5 percent of the theoretical curves. 
Thus, it is concluded that for volume ratios less than 0.2, the response of a straight 
probe can be predicted to within 5 percent by the analysis of reference 4. 

It is also of interest to determine how closely the quarter wavelength equation pre-- 
diets the resonant frequency. The quarter wavelength equation with the 8D/37T length 
correction is 


f = 



where 

f resonant frequency 

c adiabatic velocity of sound 

L length of tube 
D inside diameter 

For short probes with zero volume ratio, this equation predicts resonant frequencies 
which are higher than the reference 4 predictions by 2 to 5 percent. This agreement is 
adequate for most probe design work. If the probe has an appreciable volume ratio, a 
correction can be obtained from figure 6. 

In this work, no attempt was made to compare the experimental results with pre- 
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dictions from other existing simplified tube -and -volume analyses. It is expected that 
many of these would also be adequate for probe design work, but would be of no advan- 
tage over the simple procedure outlined in the preceding paragraph . 

From figure 6, it is apparent that a volume at the transducer end of the probe will 
decrease the resonant frequency of a straight probe. To minimize this decrease, a 
volume ratio of less than 0. 1 should be considered as a reasonable goal for a probe de- 
sign. This would result in a resonant frequency decrease of less than 10 percent from 
the ideal case. It is obvious that a shorter tube length could be used to increase the re- 
sonant frequency so long as the L/D ratio is greater than 1. However, the mounting of 
a transducer in the probe will almost always result in some volume at the end of the 
probe. Thus, as the probe length is decreased, the volume ratio will increase, thereby 
reducing, to some extent, the gain in frequency response achieved by a shorter probe. 

The frequency response of the probe with the pressure transducer mounted at right 
angles to the probe centerline is shown in figure 7. The analytical response from refer- 
ence 4 with an end correction is also shown. A comparison of figure 7 with the curves in 
figure 5 indicates that the basic response of the probe has not been altered by the right 
angle transducer mount. 


The Effects of Average Pressure Tubing on Probe Response 

The response of probes with two different sizes of average pressure tubing is shown 
in figure 8. The basic probe was the A1 -0.027; the average pressure tube was 61 centi- 
meters long and was closed at the end. One average pressure tube had a 0. 071- 
centimeter inside diameter; the other had a 0. 150 -centimeter inside diameter. The re- 
sponse of the probe with the smaller tube is negligibly different in the usable frequency 
range from the probe without the average pressure tube (see fig. 5 for comparison). The 
only apparent difference is a slight reduction in the peak amplitude ratio. The effect of 
the larger diameter tube is more pronounced. The peak amplitude ratio is considerably 
lower and the slope of the phase angle curve is reduced in the region of the resonant fre- 
quency. These changes are similar to the effects of added damping in a linear second - 
order system. One other effect is present: the amplitude ratio plot for this probe has 
small but measurable cyclic variations as a function of frequency. These variations are 
caused by standing waves in the average pressure tube interacting with the probe. 

Other tests were run with average pressure tubing in which the end of the tube was 
open to the atmosphere or terminated in a volume of 1000 cubic centimeters. The re- 
sults were similar to the closed tube results with the exception that the phase of the 
variations in amplitude ratio was reversed. Where a maximum amplitude ratio occurred 


6 



for a closed tube, a minimum amplitude ratio occurred for the open and volume ter- 
minated tubes. 

The cyclic variation in amplitude ratio for a probe with average pressure tubing can 
be predicted from the Bergh and Tijdeman equations of reference 5. This is shown in 
figure 9 where the computed response, using the subroutines in the appendix, is compared 
with the experimental data for an Al-0. 108 probe which includes an average pressure 
tube 61 centimeters long with a 0. 150-centimeter inside diameter, closed at its end. 
Figure 9(a) shows the response through the first resonance, and figure 9(b) shows an ex- 
panded plot of the response below 1 kilohertz. For the analytical curve, it was assumed 
that the average pressure tubing was connected directly into the cavity volume. The 
agreement between the predicted response and the measured data is good. Examination 
of figure 9(a) will show that successive peaks or valleys in the amplitude ratio curve are 
separated by about 270 hertz, which is roughly equal to the half -wavelength frequency of 
a 61-centimeter tube in room -temperature air. This supports the conclusion that the 
variations in the amplitude ratio of the probe are caused by standing waves in the average 
pressure tube. 

These cyclic variations in amplitude ratio represent a source of error in transient 
measurements. It is therefore of interest to determine how to limit the magnitude of 
these variations to within acceptable bounds. This might be done by reducing the diam- 
eter of the average pressure tube to the extent that viscous losses prevent the buildup of 
standing waves. The requirement for this is that D < \f2pJpTd (ref. 2, p. 229) where 
D is the inside diameter of the tube, jx is the dynamic viscosity, p is the density, and 


f is the frequency (in Hz). This approach leads to tube diameters of less than 0. 025 
centimeter, small enough to be prone to plugging. Such diameters are apparently 
smaller than necessary based on the results shown in figure 8. This is probably because 
this approach does not take into account that the pressure waves in the averaging tube 
must interact with the probe impedance to produce measurable variations. The probe 
impedance is relatively low because of large diameter of the probe. 

An alternative approach is to determine a relation between the magnitude of the 
variations in amplitude ratio and the geometry of the probe and average pressure tube. 
This was done by calculating the response of a number of different geometries and em- 
pirically finding a geometric parameter that correlates with the magnitude of the cyclic 
deviations in amplitude ratio. The geometries covered in this calculation are shown in 
table n. The calculations were done only for air at atmospheric pressure and room 
temperature. 

Figure 10 shows the results of these calculations. The average value of the peak- 
to-peak cyclic deviations in amplitude ratio over the usable frequency range (in this 


case 

eter. 


up to 600 Hz) for each geometry is plotted against an empirical geometry param- 
The parameter found to give a reasonable correlation is Dp(L^ + 13)/D^, where 
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Dp is the inside diameter of the probe, is the length, and is the inside diam- 
eter of the average pressure tube (all dimensions are in centimeters) . The line drawn 
to the right of all the data points represents a conservative estimate of the geometry 
parameter necessary to insure that the cyclic deviations in the amplitude ratio are less 
than a given value. For example, to insure that these deviations are less than 2 percent 
(i. e. , a ±1 percent error in amplitude ratio), a geometry parameter of at least 2000 is 
necessary. Thus, an Al-0 probe with an average pressure tube diameter of 0. 150 centi- 
meter would require a tube at least 125 centimeters long, but, if the tube diameter were 
0.071 centimeter, the minimum length required would be 18 centimeters. 

Also plotted in figure 10 are two measured data points for probes Al-0. 027 and 
Al-0. 108. The agreement between these measured data points and the calculations in- 
dicate that the results of figure 10 are applicable to probes with volume ratios up to at 
least 0.1. 


Effects of Entrance Baffles on Probe Response 

Test data for the baffled probe configurations are shown in figure 11, and the meas- 
ured resonant frequencies are listed in figure 3. For comparison purposes the measured 
resonant frequency of an Al-0. 108 probe is also listed in figure 3. (The volume ratio 
for the baffled probe configurations is approximately 0. 13.) These data indicate that, if 
the entrance annulus and hole areas are larger than the probe tube area, this baffle con- 
figuration does not degrade the frequency response of the probe. 

The slight increase in resonant frequency with increasing baffle area ratios is in 
qualitative agreement with the results of reference 5. (See figs. 19 and 20 of that ref.) 
The reference 5 data show that for a long tube -and -volume pressure transmission sys- 
tem, the resonant frequency can be nearly doubled by replacing half of the tubing length 
at the entrance end with a tube whose diameter is 1. 5 times the diameter of the original 
tube. This would seem to imply that similar increases in resonant frequency might be 
obtainable for baffled probes of the type tested here by optimum choice of baffle area 
ratios and relative length. Preliminary tests of this concept indicate that moderate in- 
creases (10 to 20 percent) in resonant frequency are possible. However, practical limits 
on baffle size are such that a large increase in resonant frequency is not possible. 


CONCLUDING REMARKS 

The frequency response of short transient total -pressure probe geometries has been 
measured. Variations in geometry intended to provide particle impact protection and 
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time-average total -pressure measuring capability have been studied to determine their 
effect on probe response. The measured response data have been compared with re- 
sponse predictions from various theoretical analyses. The results of this work are pre- 
sented for the purpose of facilitating the design of practical transient total pressure 
probes for turbojet engine research. 

The analyses of Iberall (ref. 4) and Bergh and Tijdeman (ref. 5) adequately calcu- 
late the amplitude ratio and phase angle versus frequency for the probes tested. The 
prediction of the resonant frequency is within 5 percent for probes with resonant fre- 
quencies up to 3 kilohertz. 

A value of 0. 1 for the ratio of transducer cavity volume to probe tube volume was 
felt to be a reasonable design goal; this volume ratio decreases the probe resonant fre- 
quency by less than 10 percent from the ideal case of zero transducer cavity volume. 

For the probes tested, the simple, well-known quarter wavelength equation was 
adequate for calculating probe resonant frequency. If used with a length correction of 
8D/3 jt where D is probe inside diameter, this equation could predict the resonant 
frequency of probes with zero volume ratio to within 10 percent. Further corrections for 
volume ratio can be made by information presented herein. 

The method of using a small diameter tube between the transducer end of the probe 
and remote steady-state pressure instrumentation in order to measure time-average 
total pressure was tested. An empirical parameter involving the probe diameter and the 
average pressure tube diameter and length was found which adequately correlates with 
interaction effects between the average pressure tube and the probe. A design for which 
the value of this parameter is greater than 2000 will insure that the error due to this 
interaction is less than ±1 percent for a probe operated at atmospheric pressure and 
room temperature. 

A baffle design for preventing particles entrained in the air stream from damaging 
the transient pressure transducer was tested. It was shown that with proper entrance 
area this baffle design does not degrade the frequency response of the probe. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 28, 1970, 

126-61. 
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APPENDIX - COMPUTER SUBROUTINE FOR CALCULATION OF THE FREQUENCY 
RESPONSE OF A MULTIPLE TUBE -VOLUME SYSTEM 


For computing purposes, a system of tubes and volumes as in figure 12 was as- 
sumed. Using the Bergh and Tijdeman (ref. 5) results as a basis, the following equations 
were used for calculation: 


D 

j-1 


V, 

cosh(^Lj) + — - (1 + sinh(<^Lj) 


where 


+ 


D M*j+l J 2< Vl W o<“j> „ 

; 1 LL (c“Sh(® )tl L j+1 ) 



-1 




1 


1 + 


y-l J 2( a i^) 
7 J 0 (“jv/Pr) 


The variables are defined as follows: 
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Engineering 

symbol 

Computer 

symbol 

Definition 

D 

D 

tube diameter 

f 

FREQ 

frequency 

J n 


Bessel function of first kind of order n 

L 

XL 

tube length 

P 


amplitude of pressure distributions 

Pr 

P 

Prandtl number (in subroutine usually equal to 0. 7) 

P s 

AM PRES 

mean pressure 

V t 


tube volume 

V v 

C 

cavity volume 

7 


specific heat ratio 


VIS 

fluid dynamic viscosity 

P s 

RHO 

mean mass density 

a 

DVDP 

dimensionless increase in cavity volume due to containment 
wall deflection 

a> = 2jrf 


angular frequency 

j 

I 

defines tube and volume 


The FORTRAN IV subroutine using the preceding equations may be referenced as 
follows: 


CALL AMPPHS(FREQ, A, PA, I,N) 

FREQ is the pressure oscillation frequency for which the amplitude ratio and phase 
angle are to be calculated and must be specified before the CALL. The program places 
the amplitude ratio P(I)/P(0) in A and the phase angle in PA. N is the total number of 
series -connected tube-volume systems. 

The system parameters must be put into the program through two common blocks 
referenced as follows: 

COMMON/ GEOMTY/D(10), XL(10), V(10), DVDP(IO) 
COMMON/THERMO/ GAMMA, P(10), RHO(IO), VIS (10) , AM PRES 
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The program accepts as many as 10 series -connected systems with diameter (D), 
length (XL), volume (V), Prandtl number (P), and volumetric displacement (DVDP) 
specified in one block. The thermodynamic properties density (RHO), viscosity (VIS), 
ambient pressure (AMPRES), and specific heat ratio (GAMMA) are specified in a second 
block. 

The program used for generating the Bessel functions is limited in use to Bessel 
functions of integer order with arguments of the form C(1 - i) where C is a real constant. 
The listing for the subroutines follow: 


SUBROUTINE AMPPHSIFREQ ,A ,PA ,1 , J) 

C 

C 

C THIS ROLTINE CALCULATES THE AMPLITUDE RATIO (A) P(I)/P(0) 

C AND PHASE ANGLE (PA) FOR A SERIES CONNECTION OF J TUBE-VOLUME 
C SYSTEMS GIVEN FREQUENCY 'FREQ*. 

C THE SYSTEM GEOMETRY IS ENTERED THROUGH COMMON BLOCK /GEOMTY / . 

C THE THERMODYNAMIC CONDITIONS ARE ENTERED THROUGH COMMON BLOCK /THERMO/. 

C 

C 

DIMENSION PRATIO(IO) 

COMMON /GEOMTY /DUO) ,XL(10 ) » VC 10) ,DVDP(10) 

COMMON/THERMO/GAMMA, P(10) ,RH0(10) ,V|S(10) .AMPRES 
COMMON /XFREQ/FREQl 
COMPLEX PRATIO, RECURS 
FREQl = FREQ 
RECUTS = 1.0 
DO 10 M= 1.10 
10 PRAT IO(M ) = 0.0 

N = J 
C 

C CALCULATE SUCCESSIVE PRESSURE RATIOS 
C 

DO 20 K=1,N 

20 CALL RATIOIPRATIO, J«-l-K, J) 

C 

C CALCULATE REQUIRED RATIO BY MULTIPLYING APPROPRIATE PRESSURE RATIOS 
C 

DO 30 K = 1,1 

30 RECURS = RECURS+PRATIO (K) 

A= CABSIRECURS) 

P A= A TAN2 ( A IMAG( RECURS) .REAL (RECURS)) 

C 

C CONVERT RADIANS TO DEGREES 
C 

PA= PA*57. 2957795 

I Ft PA .GT.0.0) PA = -360.4-PA 

RETURN 

END 
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oooooo** o s inj w oj*— n o o o n 


$ I B FTC XR AT 10 OECK 


SUBROUTINE RA TI 01 PRATI 0, J , K) 


THIS PROGRAM CALCULATES THE PRESSURE RATIO PIJJ/PIJ-l). 


DIMENSION PRAT 10 1 10) , VRATI 0 ( I 0) t CORRCTUO) 

COMMON /GEOMTY/D( 10) , XL C 10 I , V { 10) ,OVDP<10) 

COMMON / THERMO /GAMMA »P ( 10) ,RH0(10) ,VIS<10) , AMPRES 
COMMON /XFREQ /FREQ 

COM PL EX PHIt JRAT, POLY , POLRAT, CCO SH ,C SI-NH, PRAT I 0,TE RM1 , T ERM 2, T ERM3, 
CTPH I l 9 TPHI2, JRATIO 
PI = 3.1415927 

VRATIOIJ) = V(J)/IPI*D(J)**2*XLlJ)/4.) 

RFREQ = 2 .*P I*FREQ 
IF( VI J ).NE.O.O)GOTO! 

CORRCT(J) = 1.0 
GO TO 3 

CORRCT(J) = 1.0 ♦ IDVDPI J)*AMPRES/V( J)) 

CONTINUE 

IFIJ.EQ.K) GO TO 5 

POLRAT = CSQR T I POL YI J ) /POLYIJ+l) ) 

RHORAT = SQR T (RHO I J*1 ) /RHOI J) ) 

VOLRAT * D(J + l)**2/([) I J)**2I 
TPHI2 = PHKJ + 1) 

JRAT = C SQR T ( 

JJ RATIO (CSQRT(CMPLX( 0. 0,-(RFREQ*D I J + l ) **2*RH0( J+ 1 ) / 

J( 4.0*VIS(J+i ))))))/ 

JJRATIO(CSGRT(CMPLX(0. 0 , - I RF RE Q*D I J ) **2*RH0( J )/ 
j( 4#o* vis ( j mm) 

TP HI 1 = PHI(J) 

TERM1 = CCOSHI TPH I 1*XL ( J) ) 

TERM2 * VRATI 0( J)*CORRCT( J)*XL( J )*T PH 1 1 *POLY I J ) *CS I NH( T PHI 1*XL ( J ) ) 
IFIJ.EQ.K ) GO TO 2 

TER M3 = (CSINHITPHll*XLIJ) ) /CSINHIT PH 12 *XL(J+1)) )* 

T VI LRAT*RHOR A T* JR A T* POL RAT * 

T ICC0SH(TPHI2*XL ( J* 1 ) ) - PRATI 0 1 J+l ) ) 

GO TO 6 

TERM3 = CMPLXI 0.0, 0.0) 

CONTINUE 

PRATIO(J) * 1 . 0 / ( TERM1 + TE RM2 +TE RM3 ) 

RETURN 

END 


IBFTC XJRATO 


DECK 


COMPLEX FUNCTION JRATIO(CARG) 


T F IS ROUTINE CALCULATES THE RATIO OF THE BESSEL FUNCTION OF THE 
FIRST KIND, 3RDER 2 TO THE BESSEL FUNCTION OF THE FIRST KIND, 
ORDER 0 WITH COMPLEX ARGUMENT. 
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COMPLEX FM,QSUM,P SUM , XP SUM ( 2 ) ,XQSUM(2) ,XJ(2) , RAT 10, XX, Y 
IV , 7 

COMPLEX J2, JO, JR ATI 0,C ARG 
COMMON /Sl/X, Y,PSUM,QSUM 
COMMON /FMNAM/FM 

DATA P 1/3. 141 5927/, M/-1/.Z/ CO. 0,1.0)/ 

IF(REALICARG) .GT.30.0) GO TO 40 
C 

C IF THE REAL PART OF THE ARGUMENT IS GREATER THAN 30 USE IS MADE OF 
C HANKELS ASYMPTOTIC EXPANSIONS (SEE NBS HANDBOOK OF MATHEMATICAL 
C FUNCTIONS AMS NO • 55) TO CALCULATE THE BESSEL FUNCTION RATIO. THE P 
C AND Q COEFFICIENTS USED IN THE EXPANSIONS ARE CALCULATED USING THE 
C SERIES SUBROUTINE. 

C 

GO TO 20 

40 X = R E AL ( CAR G ) 

Y » - X 
DO 2 J =1, 2 
M' * M * J 
XI = FLOAT(M) 

FM * C MPL X ( XI ,0.0) 

CALL SERIES 
XPSUM( J) = PSUM 

2 XQ SIX ( J) x QSUM 
DO 3 K - 1, 2 

A = X ♦ ( 0.75 - FLOAT(K) )*PI 

XX = XPSUM(K) * COSCA) - XQSUM(K) * SIN(A) 

YY = (XPSUM(K) * SIN(A) + XQSUM(K) * COS ( A» » * Z 

3 X J ( K ) = XX ♦ YY 

RATIO = X J ( 2 ) / X J C 1 ) 

JRA7I0 = RATIO 
M * - 1 
GO TO 30 

20 CALL Z EEB E S( J 0 , 0 , C ARG ) 

CALL ZEEB E S< J 2 , 2 , CARG ) 

JRATIO = J2/J0 
30 CONTINUE 
RETURN 
END 


$ieFTC SERIE DECK 

SUBROUTINE SERIES 
C 
C 

C THIS ROLTINE CALCULATES THE VALUE OF THE SERIES NEEDED FOR 
C HANKELS ASYMPTOTIC EXPANSION 
C 

C IN CALCULATING P(J) ANO Q ( J > NUMBERS WITH EXPONENTS LESS THAN -38 
C MAY OCCUR. THIS PROGRAM USES A BUILT-IN ROUTINE WHICH SETS P( J ) 

C AND/OR 0(J) EQUAL TO ZERO WHENEVER THESE SMALL NUMBERS OCCUR. 

C 

C 

c 
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COMPLEX Z ♦ MU*PSUM,QSUM t FM f PC50I *Q(50) 

COMMON /SI /Xf Y t P SUM* QSUM 
COMMON /FMNAM/FM 
P SUM = <0.,0. ) 

QSUM = (G .,0. ) 

MU =4.*FM**2 
Z = C M PLX ( X* Y ) 

00 100 J= 1 * 50 
N = J- 1 
FN =N 

IF(N .ME. 0) GO TO 30 
20 PC J ) =f i. f 0 . ) 

Q( J ) = (MU-1. » /C8.*Z) 

GO TO SO 

30 FN 2 = <.*FN 

FN 4 = *.*FN2 
FN21 = FN 2 + 1. 

FN A 1 = 2.*FN21 

IFCCA8SCP ( J- I ) ) .EQ.O. ) GO TO 50 

PC J) = -( (MU-(FN4-3.)**2) /( 64. * (FN2-1 . I *FN2H *<< MU -< FN4-1. }**2)/Z 
$** 2 )*P ( J- 1 ) 

IF (J .GT. 2 .AND. CA8SCPCJ)) .GT. CABSCPCJ-IH ) P(JJ = (O. f 0.) 
GO TO 55 

50 PCJ ) = (O.tO.J 

55 IF( CAS SCQ ( J- 1 )) .EQ.O. ) GO TO 80 

QCJ) * -( (MU-CFN41-3. ) **2 ) / ( 64. *FN2 *F N21 ) ) *( (MU-(FN4l-l.)**2)/Z 
$**2 )*Q ( J- 1 ) 

IF (J .GT. 2 .AND. CABSCQCJ)) .GT. CABSCQCJ-IM ) OCJ) = (G.,0.) 
GO TO 50 

80 Q ( J ) - CO.tO. ) 

90 P SUM = PC J )+P SUM 
QSUM = QC J )+Q SUM 

IFCCA8SCP (J) ) .EQ.O. .AND. CABSCQCJH .EQ. 0.) GO TO 110 
100 CONTINUE 
110 CONTINUE 
RETURN 
ENO 


$ I BF TC BESSNN DECK 

SUBROUTINE ZE EBE S ( J ,N , Z ) 


C 

C THIS ROUTINE CALCULATES THE VALUE OF THE BESSEL FUNCTION OF THE 
C FIRST KIND OF INTEGER ORDER. THE ARGUMENT IS OF THE FORM 
C C( 1 - SQRTf-lH WHERE C IS A REAL CONSTANT LESS THAN 30 . 

C 

C J. ..BESSEL FUNCTION VALUE 

C N... .BESSEL FUNCTION ORDER 

C Z... BESSEL FUNCTION ARGUMENT 

C 

c 

COMPLEX J f Z 

DOUBLE PRECISION DF N , T »R ( 2 ) ,C ,D K ,D , B 
D = REALC Z) 



o o o o 


C = 0 * 1 .414213562373095D0 
CR = 0 

DFN = FLO A T ( N ) 

IF(N.NE.2> GO TO 1 
T = 5.D-1 
GO TO 2 

1 T = 1.00 

2 R ( 1 ) = T 

R ( 2 > = 0.00 
C = <C**2.) / 4.00 
OK = 0 .00 
I = 1 
M = 0 

5 DK = OK + 1.0 0 

T = T * C / ( (DK)*(DK«-DFN) ) 

I ■ I *1 

IF (T.LT.DABS(R(I)*l.D-8)) GO TO 7 
M = M ♦ 1 
B = 1.00 

IHM.GT.l ) B = -1.00 
IF(M.EQ.3) M = -1 
R(I) = R ( l ) *■ B * T 
I F( I .EQ.2 ) I = 0 
GO TO 5 

7 RR = R ( 1 ) 

R I = R(2> 

IF( N .NE.2 ) GO TO 8 
R I = (-1.) * R ( 1 ) * C 
RR =R( 2) * C 

8 J = CMPLX ( RR »R I > 

RETURN 

END 


iieFTC SUSP H I DECK 

COMPLEX FUNCTION PHIIJ) 

THIS PROGRAM CALCULATES A FUNCTION PHI WHICH IS AN 
ATTENUATION PARAMETER. 

COMMON /GEOMTY/DI 10) ,XL (10 ) , V< 101 ,DVDP(10) 

COM MON /THERMO /GAMMA ,P ( 1 0) , RHO ( 10 ) , VI S ( 10) , AMPRES 
COMMON /XFRE3 /FREQ 
COMPLEX POLY, PHI ,JRATI0 

SN =6.2831 854*FREQ*D ( J) **2*RH0< J) / VI S ( J) 

PHI =6.2831 854* FREQ* SQR T (RHO ( J) / ( GAMMA*AMPRES ) )* 

P CSQRT( l./JRATI0(CSQRT(CMPLX(0.0,-(SN/4.0)»‘ » 1* 

P CSQRT( GAMMA/P0LY( J) ) 

RETURN 
END 
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noon 


$ I BF TC SUBPLY DECK 


COMPLEX FUNCTION POLY(J) 

THIS ROUTINE CALCULATES A POLYTROPIC COEFFICIENT 
I LABEL ED *N 1 IN BERGH AND TIDJEMAN) 

COMMON /GEOMTY/D ( 10 ) ,XL ( 1 0 ) ♦ V ( 1 0) f DVDP(10) 

COMMON /THERMO /GAMMA ,P< 10) .RHO(IO) ,VIS (10) , AMPRES 
COMMON /XFR EO /FREQ 
COMPLEX JRATIOtCARGfPOLY 

SN = 6 ,28 31854*FREQ*3 ( J)**2*RH0( J) /VI S (J) 

ARG = P( J >*SN /4 * 0 

CAR G- CSOR T(CMPLX ( 0.0 f-ARG) ) 

POL Y= 1 • 0/ ( l.OH (GAMMA- 1.0) /G A MM A ) * JR AT I 0(C ARG ) ) 

RETURN 

END 


$ I BFTC FCCOSH DECK 

COMPLEX FUNCTION CCOSH(ARG) 

COMPLEX ARGtCCOSH 

CCOSH = ( CE XP ( ARG ) +CE X P ( - ARG ) )/2*0 

RETURN 

END 


1 1 BF TC FCSINH DECK 

COMPLEX FUNCTION CSINH(ARG) 

COMPLEX ARGt C S I NH 

C S IN H= (CEXP<ARG)-CEXP(“ARG) ) /2.0 

RETURN 

END 
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TABLE I. - PROPERTIES OF AIR 


USED IN COMPUTATIONS (REF. 9) 


Temperature, K 

295 

Pressure, N 

9.92 

Density, g/cm 3 

1. 172xl0' 3 

Viscosity, g/(sec)(cm) 

1. 824x10 " 4 

Specific heat ratio 

1.4 


TABLE H. - PROBE GEOMETRIES USED IN 
CALCULATING EFFECTS OF AVERAGE 
PRESSURE TUBES ON PROBE RESPONSE 


[Probe length, 2. 54 cm; probe volume ratio, 
0; average pressure tubing lengths, 25, 
102, and 203 cm.] 


Probe 

diameter, 

cm 

Average pressure tube diameter, cm 

0. 157 

0.051 

0,076 

0. 101 

.317 

.076 

,127 

. 177 

.476 

. 101 

.152 

.203 


/-Average total pressure tap 



r Pressure transducer 


Probe A 1, 0. 325 dlam 
Probe A2, 0. 160diam 


- Cavity volume 
0. 558 diam 


Figure 1. - Basic (straight) probe. (Dimensions are in centimeters.) 











Figure Z - Perpendicularly mounted transducer probe. 
(Dimensions are in centimeters. ) 



Probe 

Dimensions 

Area ratios 

Measured resonant 
frequency, 

Hz 



Annulus area to 
tube area ratio 

Hole area to tube 
area ratio 

B1 

0. 235 


0.87 

0.92 

2620 

B2 

.208 


1.30 

.92 

2770 

B3 

.177 


1.80 

.92 

2860 

B4 

.177 

0.381 

1.80 

1.80 

3030 

Al-0. 10 

— 


— 

— 

2790 


Figure 3. - Baffle probe design data. {Dimensions are in centimeters. ) 





Figure 4. - Pressure generator. 
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24 — Analytical curve with 

end correction (ref. 4)-^ 




Q> 
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Q_ 


-90- 

0 ~ 
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Frequency, Hz 

Figure 5. - Frequency response of probe Al-Ol 0 21 . 



Volume ratio, cavity volumeftube volume 


Figure 6. - Percent reduction in resonant frequency as 
function of volume ratio. Analytical curve is the 
average values computed for A1 and A2 type of probes. 
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Figure 7. - Frequency response of perpendicularly mounted trans- 
ducer probe. 



Figure 8. - Frequency response of two Al-0. 027 probes which include 
61 centimeters of average pressure tubing. 
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Figure 10. - Average cyclic deviation of amplitude 
ratio as function of geometry parameter for air 
at atmospheric pressure and room temperature. 




Figure 13. - Frequency response of baffled probe configuration. 



Figure 12. -Tube-volume system for the Bergh and Tijdemen equation. 
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